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Abstract Numerical and experimental realizations of quantum control are 
closely connected to the properties of the mapping from the control to the 
unitary propagator |14j . For bilinear quantum control problems, no general 
results are available to fully determine when this mapping is singular or not. 
In this paper we give sufficient conditions, in terms of elements of the evo- 
lution semigroup, for a trajectory to be non-singular. We identify two lists 
of "way-points" that, when reached, ensure the non-singularity of the control 
trajectory. It is found that under appropriate hypotheses one of those lists 
docs not depend on the values of the coupling operator matrix. 
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1 Introduction 

Manipulating the evolution of physical systems at the quantum level has been a 
longstanding goal from the very beginnings of laser technology. Investigations 

T.-S. Ho 

Department of Chemistry, Princeton University, Princeton, NJ 08544, USA 
E-mail: tsho@princcton.cdu 

H. Rabitz 

Department of Chemistry, Princeton University, Princeton, NJ 08544, USA 
E-mail: hrabitz@princeton.edu 

G. Turinici 

CEREMADE, Universite Paris Dauphine, Place du Marechal de Lattre de Tassigny, 75016 
Paris, France 

E-mail: Gabricl.Turinici@dauphinc.fr 



2 



Tak-San Ho et al. 



in this area accelerated after the introductions of (optimal) control theory 
tools pr, which greatly contributed to the first positive experimental results, 
see [3lil3..16,..4,..5 ,l 0lfT2] and references herein. 

Consider a quantum system in the presence of a control field e{t) G M, i > 
0. In the electric dipole approximation and in the density matrix formulation, 
the underlying time-dependent density matrix p(t) satisfies the Liouville-von 
Neumann equation 

z^p(t) = [i?o-e(%,p(t)] (1) 
p{t = 0)=po 

where Hq is the ficld-frcc Hamiltonian (including the potential) and fi is the 
dipole moment of the system. In the analysis, wc will suppose that the system 
contains N levels, thus, Hq, ji and p{t) are all iV x Hcrmitian matrices. In 
the following, for simplicity, we further assume that Hq and p arc real and 
symmetric matrices. Moreover, to avoid trivial settings, the dipole moment 
operator p is assumed to have zero trace, i.e., 

Tr{p) = 0. (2) 

In the following we adopt the notation isu{N) to denote the set of all complex 
NxN Hermitian matrices with zero trace. Recall that p{t) = U {t, 0)p{0)U* {t, 0), 
where the unitary propagator U{t, 0) G U{N), t > is the solution of the time- 
dependent Schodinger equation 

i^Uit,0)= (^Ho-e{t)pyit,0) (3) 
U{t = 0,0) = In- 

The control goal can be expressed in terms of a (self-adjoint) operator O in 
that the corresponding functional to be optimized (maximized) is (0){T) = 

Tt(^p(T)0^ . Accordingly, a maximum value of (0)(T) indicates that the con- 
trol is of optimal quality and vice-versa. 

At the heart of a good understanding of any quantum control problem, in- 
cluding the efficiency, accuracy, and stability of quantum control experiments 
as well as simulations, lies essential features of the control landscape, depict- 
ing the functional dependence of the observable on the control field, i.e. the 
mapping: e{t) {0){T). The slopes (gradient) and the curvatures (Hessian) 
are key characteristics of the control landscape. In particular, the topologies 
at and around the critical points, in which the slopes vanish, can shed im- 
portant light on the questions, for example, of what makes quantum control 
experiments (simulations) apparently "easy to perform" and why quantum 
control beats the "curse of dimensionality" (i.e., overcomes the anticipated 
exponentially growing effort required when searching over increasing numbers 
of control variables) [5]. A thorough study of the underlying control land- 
scape should give insight into the classes of future feasible quantum control 
experiments, especially those involving complex molecules and materials. 
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2 Motivation: landscape analysis and beyond 

The search procedures optimize the control quality (0){T) with respect to 
variations in the control e(t). The mapping e(t) ^ {0)(T) is the composition 
of two maps: a control map e{t) i-^ U{T, 0) and a "kinematic" map U{T, 0) i— ?• 
(0)(T). A key issue is whether there exist critical points that are suboptimal 
solutions. To analyse the critical points it is expedient to consider these two 
maps separately. 

Previous analyses |14| have investigated the mapping from the propagator 
to the observable i.e. U{T,0) n> (0)(T); the critical points correspond to the 
equation 

S{0}{T) 



6U{T, 0) 



0. (4) 



Assuming controllability hypothesis pi[T5] , this equation is satisfied if and only 
if 

[U*{T,0)OU{T,0),pm^0, (5) 

i.e, we obtain a necessary and sufficient condition for the kinematic critical 
points [7] . The optimization " lanscape" , when looked upon in terms of this so 
called "kinematic" mapping, is therefore exempt of suboptimal critical points. 

To extend the landscape analysis beyond the kinematic setting, i.e. to in- 
vestigate the dependence on the control itself, we need to analyse the solutions 
of the following critical point condition [S] 

The solutions of this equation are among the critical points of the control map 
e(t) i-^ C/(T, 0); it was shown in |17l that this mapping is nonsingular when all 
matrix entries of the dipole moment 

m = u*it,o)^iUit,o), (7) 

as functions of t £ [OjT'], are linearly independent (to the extent possible), or 
equivalently if the dipole moment matrices /i(i)'s span the space of all zero- 
traced Hermitian matrices isu{N). 

In previous works, it has always been assumed that "full controllability" 
implies "fully linear independence" of the dipole moment matrix entries. An 
assertion of the linear independence of dipole moment matrix entries is also 
essential for understanding the landscape of unitary transformation control 
problem [8]. In this paper, we present a rigorous analysis of this assertion 
and identify properties of the trajectory, in terms of the propagators, that 
are sufficient to ensure the linear independence of entries of fi{t). In addition, 
the analysis of which controls correspond to a non-singular mapping e{t) i— 
U{T, 0) or €{t) t-^ (0)(r) is also relevant, e.g., for stabilization purposes, error 
cancellation, maintaining coherence, etc. 
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3 Dipole dependent way-points 

Denote by L the Lie algebra generated by — i-ffo and — i/i as a sub-algebra 
of u{N). We will suppose that the system is density matrix controllable, or 
equivalently 

L = su(7V) or L = it(7V). (8) 

Note that when the system is density matrix controllable, then it is also wave 
function controllable, cf [BE]- 

Theorem 1 Under assumptions Q and ^ there exists a set {Ui, U2N'^-2n} 
of unitary propagators (that we will call the set of "way-points") such that if 
U{t,0) visits all of them i.e. U{tk,0) = Uk for some times tk, fc = 1, 2iV2 _ 2N 
the components fij{t) of the matrix fi{t) in ^ are linearly independent as 
functions of time. 

Consequently, under the hypotheses 1^ and ^ for T sufficiently large 
a control field e{t), t 6 [0,T] exists such that functions fij{t) are linearly 
independent (or equivalently .such that fi{t) spans the set of all zero-traced 
Hermitian matrices). 

Remark 1 Note that here there is no information required on the structure 
of the matrices Hq or jjL. 

Proof Since the Hermitian matrix /i has zero trace it cannot be a constant 
and thus it has at least two different eigenvalues Ai and A2. 

Let us begin by some general conventions: for any i ^ j and any matrix 
M we denote by M^ij^ its 2 x 2 sub-matrix formed by the i-th and j-th rows 
and columns: 



M. 



<h3> 



(9) 



Moreover, when we alter the N x N identity matrix by putting in the i-th and 
j-th rows and columns a given 2x2 matrix D, we will denote the resultant 
N X N matrix by D<'J> : 



]J<i,3> 









J 
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We choose Uk in the following way: let Z = 1, . 
= l,...,iV with i < j. 

'Al 
A2 



U41+1 is such that (fi(t4i+i) 



<i,j> 



(10) 



^''^y index the couples 
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C^4;+2 = U41+1 f 20) ' ^^'^'^ (^KUi+2)j = { " I and all 



<l,J> 



other entries are identical to those of //(t4;_|_i). 
- C^4;+3 5C/4/+1 ( 2^ ) ■ Then 



^ 1 / Ai + A2 Ai - A2 
<i,j> 2 y Ai — A2 Ai + A2 



(11) 



and all other entries are identical to those of fllt^i^i). 

<i,J> 



Uii+i = \Uii+i ( 2 j ) ■ Then 



A 1 / Ai +A2 (A2 - Ai)A , . 

^'^'^'+*))<..> = H(Ai-A.)^ A,+A2 j (^2) 

and all other entries are identical to those of /i(t4/4-i). 

Note that all /i(tfc) = U^fiUk are Hermitian (and of zero trace). Let us take 
a zero trace Hermitian matrix Z such that 

=0, Vfc = l,...,4iV. (13) 

where we used the canonical scalar product of Hermitian matrices (A, B) = 
Tr{A* B) = Tr{AB). We recall also the definition of the norm induced by this 
scalar product ||^|| = ^ {A, A). 

We denote (^^^ = Z<jj>; since {Z,ft{tk)) for fc = 4/ + 1, 4/ + 4 

it follows that {Z^ijy, (/*(ifc))<jj">) is the same for fc = 4^ + 1, Al + 4. The 
consequence is that A1X+A22: = A2.i; + Aiz which together with Ai ^ A2 implies 
z. Moreover we also have Ai2: + A2a: = ^ ^2(Ai + A2)x+ (Ai — A2)(y + y*)^ 

thus y + y* =0. Additionally, we have Xix + X2X ~ i^2(Ai + -^2)2; + (Ai — 

^2)i{y* — y)j , thus y ~ y*. From y + y* =0 and y = y* we infer that y = 
and obtain Z^ij^ = XI2 {I2 is the 2x2 unit matrix). Since this is true for 
arbitrary i < j we obtain Z = xIn- But, since Tr{Z) = it follows Z = 0, 
q.e.d. 

For the second part of the conclusion we use the fact that under the hy- 
potheses ^ and ^ the system is controllable thus for T sufficiently large 
a control field e(i), t £ [0,T] exists such that functions fij{t) arc linearly 
independent. 



X 



4 Dipole independent way-points 

The result of the previous section can be interpreted as follows: as long as 
the control field ensures that the propagator will visit some specific unitary 
transformations, then the required linear independence of the time-dependent 
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elements of fi{t) is satisfied. Note that this set of unitary transformations (i.e., 
propagators) expUcitly depends on the dipole moment. In the laboratory, pro- 
cedures have been designed that are able to control the quantum evolution 
even in the absence of precise information on the dipole [TT]. It is therefore 
possible in principle to experimentally find a control passing through a spec- 
ified list of propagators when this list does not depend on the dipole entries. 
A universal procedure can therefore be implemented that will ensure the non- 
singularity of the mapping e{t) i— > U{T, 0) even when the dipole is unknown. 
Beyond the question of the landscape analysis, such a procedure can be useful 
in additional circumstances when the mapping is required to be non-singular, 
e.g. for stabilization purposes, etc. We will therefore investigate the following 
circumstance: suppose that the system is controllable and thus one can exper- 
imentally implement arbitrary propagators U G SU{N). Can we find a list of 
propagators (as in the previous section) which are independent of the precise 
values of entries of ? 

Theorem 2 Under assumptions ^ and ^ there exists a set of "way-points" 
W C SU{N) independent of fi such that if U{t, 0) visits all propagators in W 
(i.e. for all [7^ 6 W there exists tk with U{tk,0) ~ Uu) the components fij{t) 
of the matrix fi{t) in (jT]) are linearly independent as functions of time. 

Proof We denote for any U G SU{N): 



Cu = {(Z,^) e {isu{N)f I \\Z\\ = IIa^II = 1, Tr{ZU*^iU) ^ 0}. (14) 



This set is open because its complement is the solution of a linear equation in 
the entries of the matrices Z and fi. 

We prove now that for any Z, /i £ isu{N) of unit norm there exists a 
U e SU{N) such that Tr{ZU*^iU) ^ 0. To see this, diagonalize Z and n: 
Z = UlDiUi and /i = U2D2U2 with Ui,U2 unitary and Di,D2 diagonal; we 
will denote by dJJ the k-th diagonal entry of . Recall that since Z, /i are 
zero-traced rf^ = for a = 1, 2. Therefore, in each set ((ifc)^i there are 

some elements that are strictly negative and some strictly positive. We will 
also suppose that the diagonalization is done in such a way that the diagonal 
elements are ordered from the lowest (that is strictly negative) to the highest, 
which is strictly positive. Then, 



Tr{ZU*fiU) = Tr{UlDiUiU*U;D2U2U) = Tr{DiV*D2V) (15) 



with V = U2UUI G SU{N). Take now U such that V is a permutation matrix 
corresponding to permutation a. Then 



If X]fc=i ^k^a{k) ~ ^ permutation a then it follows in particular that 

SfcLi ^fe^fe = '^i'^w-fc' ^^"^^ ^^^^ member is always superior to the 



N 




(16) 



fe=i 
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second with the equality implying that one of the vectors is constant which 
leads to a contradiction with the zero-trace hypothesis. 

We therefore obtain that ^u£SU(n)Cu is a covering of the compact set 
{{Z^^) £ {isu{N)Y I \\Z\\ ~ W^W = 1} and thus one can extract a finite 
covering i.e. a finite set W such that 

Uc/ew Cu - {{Z, m) e {isu{N) f \ \\Z\\ = = 1}. (17) 

This means that for any Z and /z in isu{N) (note that there is no need to 
impose the norm condition) there exists U € W such that Tr{ZU* jJU) ^ 0. 
Expressed otherwise, for any /i 6 isu{N) there does not exist Z G isu{N) such 
that Tr{ZU*nU) = for any t/ e W. This implies that | [/ G W} 

form an independent set of vectors in isu{N), q.c.d. □ 

Remark 2 Theorem 2 only assures that such a set W exists but does not give 
yet precise information on its size (which can be large) nor does it provide a 
constructive approach to identify in practice W. The following result addresses 
these questions upon imposing an additional constraint on the dipole fi. 

Theorem 3 Under assumptions ^ and (|5]) and 

^ 0, Vi ^ j. (18) 

there exists a set of "way-points" § C SU{N) such that if U{t,0) visits all 
propagators in § (i.e. for all [/fc G § there exists tk with U(tk,0) = Uk) the 
components fij{t) of the matrix jlit) in Q are linearly independent as func- 
tions of time. Moreover the set § can be chosen to be the same for all coupling 
operators /i satisfying the hypotheses 

Proof Suppose that the entries of jl{t) are not linearly independent i.e. a 
matrix Z , Tr{Z) = exists such that Z is orthogonal to all matrices U*^U, 
yU G S, i.e., Tr{ZU*pLU) = VJ7 G S. 

Let us compute U* : for ^ = i ^"""^ ) (here is a 2 x 2 matrix, ^22 



a N -2 X N -2 matrix, etc.) and U = W 



— w<i,2> _ 



W 

In-2 



U^.U^r^^lf'^^^^A. (19) 

A typical example of W is | ^^'^^f^.s I . The entries of W*niiW will 

^ \^sm((/)) - cos(0)y ^ 

contain terms of second order in cos(0) and s'm{(f>): cos^(0) = cos(20)+i 
sin2(0) = i^£^(M)^ sin(</.) cos(0) = SIBMI- the entries of W^>i2 and H21W 
will contain only first order terms in cos((/)) and sin((/)) while /i22 docs not 
depend on (p at all. 

We recall that for any 6* G M the functions 1, cos(6'), sin(6'), cos(26'), sin(26') 
arc linearly independent, which will be used in the following fashion 
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Lemma 1 There exist five real values of 9, for instance ~ {0, 7r/3, 7r/2, tt, Stt/S} 
such that 

OL\ + OL2 cos{9) + Q!3 sin((?) 

+a4 cos(26l) + as sin(26') = Ve* e (20) 

implies 

ai = ... = a5 = 0. (21) 

We now introduce into the set § the matrices 

_/ cos(0)+*sin(0)\<''^> 

^''■'-■^ ~ \cos{9)-ism{e) j ' ^ ' 

for any i 7^ j and 9^0. 

Since Tr{ZUg ^.jfiUe.i-j) = V6' S 6* by Lemma[l] the coefficients of sin 20 
and cos 20 in Tr{ZU* ^U) = will be zero. The coefficient of sin 20 turns out 
to be fiij{Zij + Zji) and that of cos 20 is iiiij{Zij — Zji). It follows from 
that Zij = Zji — 0. 

Let us add to the set § the matrices 

_ /cos(0) sin(0) 
'^'■'•^ ~ \ sin(0) - cos(0) ' ' ^^'^^ 



for any i 1, iV-1, j = i + 1 and £ 6>. We use again Tr{ZV^^. -[iVe^i:f) = 
\/9 € and identify the coefficients of sin 20 and cos 20 which arc respectively 
Hij{Zu - Zjj) + ' " 2''" {Zij + Zj.i) and —^ij(Zij + Zji) + ^" 2' " (^ii - ^jj)- 
We obtain = Zjj. Since this is true for any i = 1, N — 1 from Tr{Z) ~ 
it follows that Z„- = Vz < iV. 



Remark 3 Although the hypothesis (|18p maj/ seem strong we believe that it is 
rather a technical requirement which could hopefully be relaxed in the future. 



5 Conclusion and discussion 

Following previous works in [T3] which study the properties of the control 
input-output map we investigate in this paper the singularity of the control to 
propagator map through the study of the linear independence of entries of the 
time-dependent coupling operator fi{t) (cf. (O). We provide several sufficient 
conditions in terms of the evolution semigroup points that are reached during 
the propagation of Eq. ([T]) . We expect that the criterion in ([T5|) can be relaxed 
to accommodate even more general coupling operators. 
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